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ABSTRACT 



An analysis of the magnetic effects of the quiet time proton belt 
has been performed with data obtained from a proton detector aboard 
Explorer XII, From the measured energy spectra, intensities, and pitch 
angle distributions of the protons with energies above 100 kev, the current 
density on a magnetic meridian plane was calculated using the theory of 
Akasofu and Chapman, modified to eliminate several simplifying assumptions. 
The algebraic sum of all the currents was 0.59 million amperes, and the 
magnetic moment of the current loop was 0.0291^. From the electric 
current distribution the magnetic perturbations on the meridian plane 
were obtained. At the magnetic equator on the earth's surface the 
proton belt produced a 9 y decrease in the field in comparison with 
the 38 y decrease reported by Akasofu, Cain and Chapman based on some 
preliminary data from the same detector. The maximum perturbation 
appeared at L = 3.6 on the magnetic equator, where the field attained 
a value of -23 y. 
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I. Introduction. 


The various motions which a group of charged particles exhibit when 
trapped in the earth's magnetic field result in electric currents. The 
analysis of these movements yields the observation that the net currents 
are parallel to lines of magnetic latitude: in some regions of space they 
flow eastward, and in other regions westward. If one assumes longitudinal 
symmetry about the earth, these motions necessarily correspond to a toroidal 
ring current encircling the earth. 

On the basis of geomagnetic evidence it has long been postulated that 
some type of ring current exists. The enhancement of the current at times 
was thought to be one cause of magnetic activity, but even during 
magnetically quiet times the current was never supposed to die away 
completely. 

Thus, since the advent of the satellite era and the discovery of the 
earth's trapped radiation, there has been considerable expectation that the 
particles which constitute the ring current would be directly observed, and 
their species, energies and spatial distributions determined. In addition, 
magneticians have sought the location of the currents by attempting to 
measure their magnetic effects with satellite borne magnetometers. While 
the trapped radiation was discovered over five years ago, the endeavor to 
observe any significant ring current has been unsuccessful until recently. 

In an initial survey of data received from a low energy proton detector 
aboard Explorer XII (launched August 16, 1961), Davis and Williamson (1963) 


noticed that the energy density of the 100 kev to 4.5 Mev protons amounted 
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to an appreciable fraction of the energy density of the earth 1 s magnetic 
field during a magnetically quiet period, and from 4 out to 8 earth radii 
the ratio of the two energy densities was almost constant* It was of obvious 
interest to determine the magnitude of the magnetic field produced by this 
population of protons. On the basis of some preliminary data furnished by 
Davis, Akasofu , Cain and Chapman (1962) calculated the magnetic disturbance 
for a model proton belt which resembled the measured proton distribution. 

The results predicted a decrease in the surface equatorial field of about 
38 y due to the quiet time belt. Davis and Williamson (1963) also observed 
a three fold increase in proton intensities at small pitch angles during a 
magnetic storm. Adopting the assumptions that (1) this increase was 
applicable to protons of all pitch angles, and (2) the average energy of 
the particles equalled the quiet time energy, Akasofu and Cain (1962a) 
calculated that an additional 80 y decrease would be observed at the 
equator of the earth, sufficient to account for the -60 y main phase of the 
aforementioned magnetic storm. 

A considerably more thorough analysis of the magnetic effects of the 
quiet time proton belt has now been performed from the data obtained from 
the proton detector aboard Explorer XII, and the results are transmitted 
in this paper. In addition to the fact that these protons produce the 
greatest magnetic disturbance of any known population of trapped particles, 
it is also of interest to study the effects of the ring current on 
phenomena related to the earth’s main magnetic field. The magnetic lines 
of force are distorted from the shapes calculated by the spherical harmonic 
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analysis of the field at the earth's surface. There is a stretching out or 
strengthening of the field at distances beyond the ring current which affects 
the location of the magnetospheric boundary. Of interest to cosmic ray 
physicists is the influence of the Stormer cut-off rigidity for solar protons 
approaching the earth near the polar regions. Finally, the quiet time 
particle distributions are established in order to acquire a basis for 
determining storm time variations. 

Thus, in this paper the proton data are also interpreted so as to 
determine the extent which the magnetic perturbations of the protons effect 
the earth's main field. 

II. Proton Lifetimes. 

The proton ring current under discussion has been described as a "quiet 
time" current, there being arguments against consideration that its 
enhancement during a magnetic storm is the cause of the main phase. If 
one assumes that the lifetimes of these particles are limited solely by 
coulomb scattering and charge exchange, the decay time of an enhanced 
current during a magnetic storm would be an order of magnitude larger than 
the one day recovery time indicated by magnetograms at ground level ( Liemohn , 
1961 ) 0 To obtain agreement between the observed decay time and calculated 
proton lifetimes, it is necessary to postulate that protons in the energy 
region 1 to 100 kev (below the threshold of the Explorer XII detector) 
constitute the storm time ring current ( Dessler , Hanson and Parker , 1961). 

In light of the fact that at magnetically quiet periods the observed proton 
intensities remain stable to within a factor of two over a long period of 
time ( Davis , Hoffman and Williamson , 1964), the electric current resulting 
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from them is considered to exist at all times and is described as a quiet 
time ring current. 

However, one must recognize the fact that an enhancement of the protons 
within the energy range of the detector might constitute the higher energy 
portion of the storm time particles. Since protons trapped at a given time 
do have finite lifetimes, there must be a source and probably an acceleration 
mechanism for replenishing the population* In view of the insufficiency of 
our knowledge concerning such processes, it is possible that the particular 
acceleration mechanism, and perhaps others, may influence the lifetimes of 
the particles. It has also been experimentally observed that during a 
magnetic storm some time variations do exist. For these reasons the storm 
time proton observations are being carefully analyzed and will be reported 
in a later paper. 

III. Electric Currents . 

Since the magnetic effects of the trapped protons are due to the 
electric currents they establish, it is first necessary to deduce the 
current distribution on a magnetic meridian plane from the experimental 
observations of the protons. For this calculation a slightly modified 
version of the theory developed by Akasofu and Chapman (1961) is used. 
Essentially the same notation has also been adopted and appears for reference 
in the Appendix* 

In Alfven f s guiding center approximation describing the motion of 
charged particles in a strong magnetic field, the motion is analyzed into 
three parts: (a) the gyration about and in a plane perpendicular to a 
magnetic field line; (b) the oscillati on essentially along a line of 
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force between two mirror points; and (c) the east -west drift motion for 
positively charged particles on a magnetic shell. All of these motions 
may contribute to the flow of an electric current, i. To calculate the 
current produced in a gas permeated by a magnetic field, the gas pressure 
and the characteristics of the magnetic field must be known at every point 
in the gas. 

We now limit the discussion to the steady state case with a proton 
distribution symmetrical about the axis of the magnetic poles. While 
the experimental data will be analyzed in the B/B g (or equatorial pitch 
angle), L coordinate system, the computations will assume a geomagnetic 
dipole field. If lc is the eastward unit vector normal to a magnetic 
meridian plane at a particular point, the volume current densities arising 
from the three motions are obtained from the following formulae (Parker . 
1957): 

a) The current produced by the gyration: 


tL ‘ k [nKJ (SK* t it-] 


where 

h = cos cp = 1 / B \ 

2 (1+3 sin^cp)^ cos-^cp / 

and 

— = 3 (1 + sin^tp)- 

R c aL coscp (1 + 3 sin^cp)^^ * 

b) The current due to the oscillations: 


-1 


( 1 ) 


( 2 ) 


(3) 


i o = enw s = °> 


(4) 
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because as many particles are passing through a unit area perpendicular 
to B towards one mirror point as are moving towards the other mirror 
point, at least over a short period of time, 
c) The current produced by the drift motion: 


A 

i = k 
D 



Eli<>Pm 

Pm nr 



( 5 ) 


The total current is then 


A 

- ik 


with 


i 



C iPn 
aBh 2 >- L 


( 6 ) 

( 7 ) 


The negative sign in equation 6 signifies that a positive term in i 
is a westward contribution to the current density, and a negative term is 
an eastward contribution. 

Let us briefly consider just what details in the motions of the 
particles give rise to the three current terms of i. The first term in- 
volving p s is due to the westward drift motion of protons which depends 

2 

upon the centrifugal force mw g /R Q associated with the motion along a 
curved field line. 

The second term also arises from the curvature of the lines of force 
The circular orbits due to the gyration are more crowded together on the 
concave side of the line than on the convex side. This is illustrated in 
Figure la, where between two adjacent lines of force there is a higher 
density of particles per unit area moving into the page than out of it. 
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For protons this results in a net eastward current at all points. 

Surprisingly, this term is larger than the term due to the drift motion 
if the pitch angle distributions are reasonably steep, as they are for 
the protons in the belt being discussed. 

The final term is due to the fact that the radiation exists in a belt: 
the intensity increases with radial distance to a certain point, and then 
decreases. The effect in producing currents is illustrated in Figure lb, 
in which the particle intensity is larger on the line of force at a greater 
distance from the earth, causing a more intense movement of charged particles 
in a direction into the page than out of it in the region between the lines 
of force where the two currents intersect. On the inner side of the belt 
the current will be eastward, and on the outer side it will be westward. 

The calculation of the current at all points on a meridian plane 
first requires the determination of the pressures p g and p n at each point. 

To simplify the equations we first consider protons of a single speed w. 

2 

The longitudinal pressure p s * = Nmw g . 

The lateral pressure p n * = ^Nmw n • 

While it would be possible to calculate these pressures at each point 
in the magnetosphere from the experimental data, considerable simplification 
is obtained if they can be expressed in terms of a few parameters determined by 
the particle distributions. This is accomplished by first fitting the data to 
a special equatorial pitch angle distribution function with which the local 
intensities and pitch angle distributions can be obtained from the equatorial 
intensity and distribution by means of a single parameter depending only upon 
L. Then the current equation is integrated over the energy spectrum 
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of the particles of interest at each L value; this integral turns out to 
be actually the total momentum flux at the equator* Hence, the local 
current at a given latitude depends only upon the two parameters 
representing the pitch angle distribution and equatorial momentum flux. 

Let F(L,s,9) denote the number density of protons at a point s on 
a magnetic shell L, whose speeds lie in dw at w, whose pitch angles lie 
in d0 at 0, and whose azimuth angle is arbitrary, since we assume 
symmetry in azimuth. Note that this definition of F(L,s,0) is slightly 
different than that of Akasofu and Chapman (1961) in that it includes a 
differential angle in azimuth. Fortunately, we can also adopt Parker's 
special pitch angle distribution ( Parker . 1957) for the functional form 
of the distribution in pitch angle, but with the conditions discussed in 
Section 4 "Proton Data" of this paper: 

v 2 

F(L,S,e) = Cq, / ^e\ sin a e. (8) 

\B i 

This distribution has very useful and curious properties. It can be 

divided into two parts: the second part, sin 0, 0, is the local pitch angle 

distribution and has a form independent of B e /B. Hence the ratio of p n to 

p remains the same at every point on an L shell. Because the form function, 

sin a 0, is independent of B e /B, the integral over all pitch angles, which 

would be evaluated to obtain the omnidirectional intensity at any value of 

B p /B, is independent of this ratio of fields. Hence, the omnidirectional 

a 

intensity varies down a line of force only as (B e /B)2". 
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The constant can be expressed in terms of the total density at 


the equator N e (L): 


N (L) = C 
e cf 


r " 

sin* 0 (2 tt sin 0) d0 
0 . 


(9) 


The integral in N (L) will arise again in the theory and can be 
e 

expressed in terms of gamma functions: 

~ tt 


VO 


( a+2 \ 

sin** 1 © d0 =y^ 2 


TW) 


( 10 ) 


Hence 


C 

a 



TT 


-3/2 



N e (L) = A(«)N e (L) 


Therefore 

F(L,s,0) = N (L) A (a) 



sin 


a 


0 . 


(ID 


( 12 ) 


Since a particle detector measures a flux rather than a density it 

would be well to express all equations in terms of fluxes. 

or 

F(L, S> 0) = I e< L ) A(q)| B ej sin* 0 (13) 

where I e (L) is the omnidirectional intensity at the magnetic equator as a 
function of L for protons of speed w. It must be noted that a is a 
function of L, and could also depend upon w. 

Inserting this expression into the formulae for the gas pressures 


- 10 - 


yields for protons with speed w: 


n 


TT 


p * = m 
r s 


a_ 

2 


(w^cos^G') x e(L) a(q/) ( ®e\ sin a 9(2 tt sin 0) d0 (14) 

\ B“ 


w 


and 


r>. rr 


P n * = E 


(JO 


w 


(w 2 sin 2 0) x e ^ L) A (a) / V\ sin a 0 (2 rr sin 0) d0 (15) 

\ B / 


Evaluating the integrals with the aid of equation 10, we have 


p n * = I mwI e (L) n (o?)^^£ j 


a_ 

2 


(16) 


where 


n (<*> = ag 


and 


P 


* 


s 



P *. 
y n 9 


(17) 


(18) 


which indicates again that the ratio of partial pressures is constant on an 


L shell. 
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The current can now be expressed in terms of I fi (L) and a- 


i* = c a 


BR 


c a+2 


aBh, 2 


+ 


*- - c 

3Pn* 

n aBh2 

3L 


L mwl (L) ( 

!e\ 

a 

2 

2 6 V 

B / 


[ H(a)/Be 

Of 

-f 

dI e (L) 

L V B 


dL 

a 

( — ^ ln ( 


\l9L ] 

\B / ‘ 

* B 1 

/ 9L J 


fr) 


a 

2 


I e (L) 


(19) 


or 


i* = 


2 3 3 

mwc f D*(o',tp)L I (L) + E*(of,cp)L ^£e(L) + F*(o'>cp)L I (L) 5a(L)"l (20) 

«- *L H J 


where H is the field strength on the equator at the earth's surface, and 
o 

or j_ / 

2+4 


D*(a',tp) = \ I !±\ _a_ ii+siruii , 

Mb/ *+3 col i 9<(> 


E*(a,cp) = 1. <y+2 / B e\ 
2 a +3 \ B / 


COS cp , 


( 21 ) 


F*(o;»cp) — E*(o;,cp) T 1 — 1_ In (— \ 

1_ (or+2)(ar+3) 2 \ B e/_ 
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Equation 20 corresponds to equation 52 of Akasofu and Chapman (1961) 
with the added term depending upon a (L). 

It is now necessary to integrate the current over the energy range 
measured by the detector. Since the scintillator has a low energy cut- 
off near 100 kev, the current density distribution arising from protons 
with energies greater than 100 kev will be calculated. This low energy 
limit must be kept clearly in mind in what follows. 

We next make another simplifying assumption: that oi is independent 
of energy, or the energy spectrum is a function only of L, not of pitch 
angle. The validity of this assumption will be discussed in Section 4, 
"Proton Data." This allows D*, E*, and F* to be independent of energy. 

The only terms remaining in the current equation 20 which depend 

upon energy are the velocity w and the equatorial omnidirectional flux 

of this velocity, I(L), which is actually the omnidirectional differential 

energy spectrum, ^L. Let 
dE 



which is a momentum flux, or twice the energy density of the protons* 
Then 


i ■= 


C r D(a,cp) lVl) + d*(L) + F(a,<p)L 3 *(L) 

aH I dL dL 

o 


(23) 
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Hence, the current at any point in space, L, cp> need be expressed only in 
terms of the pitch angle parameter a and the momentum profile f at the 
equator. 

IV. Proton Data. 

The data relating to the low energy protons from which the current 
distributions have been calculated were acquired from one orbit of the 
satellite, pass 10A (outbound) and 10B (inbound) on August 27 and 28, 

1961. At this time the major axis of the orbit was lying in the 
meridian plane which included the sun (see Figure 2). Since the heart 
of the proton belt (L = 3.5) lies at a relatively low altitude in com- 
parison with apogee (L = 13) much of the data were collected near the 
twilight regions, although the outer edge of the belt was traversed on 
the sunlit side of the earth. 

The planetary magnetic three -hour -range index Kp, plotted in Figure 3, 
indicated that no major disturbances were in progress. In fact the first 
appreciable magnetic disturbance during the lifetime of Explorer XII did 
not occur until pass 12. 

This particular pass was chosen for thorough analysis because the 
orbit lay closer to the geomagnetic equator in the heart of the belt than 
any other orbit. Also the telemetry coverage was complete. When this 
study was initiated it was expected that a very typical proton dis- 
tribution would be observed during the pass 10, since Davis and Williamson 
(1963) found no temporal variations larger than 307« during nineteen passes 
of Explorer VII distributed before and after this pass, but not including 
it. However, later comparisons with this original data indicate that the 
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proton intensities were enhanced by as much as a factor of two at L values 
smaller than the maximum intensity at L of 3,5 during pass 10. The 
maximum was also somewhat sharper than observed on the average pass with 
a more rapid decrease in intensity with increasing L beyond the maximum. 

This means that the belt was lying at slightly smaller L values than the 
average belt. 

The details of the low energy proton detector have been described by 
Davis and Williamson (1963). Let it suffice here to produce a table containing 
the spectra employed in calculating the momentum integral jf (Table 1), as 
well as some examples of pitch angle distributions. In general the energy 
spectrum was described as a doubly sloped exponential spectrum, illustrated 
in Figure 4 in the integral form which was actually measured by the 
detector. At L values below 2.2 the data points were best fit with power 
law spectra, and up to L of 3 a steep power law spectrum was included for 
the higher energy protons, as suggested by the data of Bame et al (1963) 
and Naugle and Kniffen (1963). The spectra listed in Table 1 were 
differential before insertion into the momentum integral. 

In the derivation of the formulae from which the current distributions 
may be computed, considerable simplification was attained by employing 
as the functional form of the pitch angle distribution, sin^ 0. The 
measured angular distributions were found to fit this form reasonably 
well provided the entire range of pitch angles, 90° to 0°, was divided 
into two regions, 90° to 30°, and 30° to 0°, each with an independently 
determined value of cv* Some sample fits of the data to this function appear 


in Figure 5. 


-15- 


Figure 5a illustrates the extrapolation to 90° of the distribution 
based on data points with pitch angles only between 30° and 56°. At 
low L values the orbit carried the satellite away from the magnetic 
equator making the measurement of protons with large pitch angles 
impossible* Figure 5b is an excellent example of the necessity of using 
the two independent values of or. An attempt at a least squares fit of 
all data from 0° to 90° produced a curve (solid line) with a large hump 
at the middle pitch angles and a very low value at large pitch angles* 

Another assumption which considerably simplified the derivation of 
the current distribution formulae involved the independence of a upon 
energy. In reality the value of a increases as a function of increasing 
energy, as illustrated by the data in Figure 6 for protons of equatorial 
pitch angles between 90° and 30° • The data in the figure were obtained 
from the output of the detector which measured the incident energy flux of 
protons rather than the particle flux* The energies at which the points 
are plotted are the effective low energy cut-offs for the integral energy 
spectrum measured at this value of L. By choosing a constant value of ot 
for a given L, the functions D*, E*, and F* (equation 21) became energy 
Independent and could be removed from the integration over energy* Hence, 
the value of a (L) selected is approximately the value derived for the 
particles having the average momentum on the given L shell* This value of 
a as a function of L is plotted in Figure 7 for each of the two pitch 
angle ranges, 90° to 30°, and 30° to 0°. 

The calculation of employed in the determination of the momentum 

integral ^ (see Figure 4 for spectral shape of the integral), required 
some special attention* This quantity, the Integral omnidirectional 
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intensity above zero energy, is equivalent to the directional intensities 
at 100 kev integrated over all pitch angles, which can be actually measured 
by the detector, multiplied by the factor exp i 100/F. Q J , However, the 
counting rates of the lowest energy channels of the detector aboard 
Explorer XII passed into saturation about one earth radius on either 
side of the maximum intensity at 3.5 , at least at large pitch angles. 

Hence, it was necessary to compute by either of two less direct methods. 

If (1) a directional intensity above an energy E , J(> E ), is measured 

o c 

at a pitch angle 0 C , (2) the pitch angle distribution characterized by the 
parameter a is known, and (3) the energy spectrum (E q ) is measured (the 
latter two parameters could be determined by the current channels of the 
detector, which did not pass into saturation), then the value of N can 
be calculated from the formulae 


Nj (> E) » J(> Ec ) 

C 0 sin® e c 



r <v e >/e„] 


(24) 


where G q is the telescope factor = 5.85x10 cm -ster., the expression 
2r?^ times the gamma functions arises from integration over the pitch 
angles, (sin** 0 )“* is actually an extrapolation of the directional Intensity 
measured at 0 C to the intensity at 90° equatorial pitch angle, and 
exp f (E C -E)/E Q j is an extrapolation in energy from the integral point 
at E c to the energy E. The quantity Nj is then the integral omnidirection 
flux above the energy E, so N would have E = 0. 
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Therefore, when the counting rate channel at 100 kev was in saturation 
either of two extrapolations could be used to acquire A higher energy 

channel (for example, the 880 kev step) would provide a non-saturated 
counting rate at a near 90° pitch angle, making (sin a 0 )”^ — 1 but 
requiring an extrapolation in energy exp [ (E c -E)/E q J , or a smaller pitch 
angle measurement from the 100 kev channel would provide a measurement of 
J(> E c ) out of saturation, allowing exp £ (E c -E)/E 0 ]| to be near one, but 
increasing the value of (sin a 0 )“^. Obviously the accuracy of either 
method depends upon how well the data can be characterized by the parameters 
O' and E o » The values of Nq^ employed in the computation of % (L) were 
acquired by utilizing both methods. Table II shows a sample set of 
intensities calculated for an L value of 3.25 for particles with pitch 
angles 30° to 90°. The agreement between the various calculations of an 
at a given L value was within + 15%. 

From the two sets of values of a (L) necessary to describe the protons 
over the entire pitch angle range from 0° to 90°, a double set of is 
calculated, since this quantity is a function of a through equation 24. 
Hence, we obtain two independent profiles of (L), one to be associated 
with protons having pitch angles 30° to 90°, the other with protons from 
0° to 30°. These two profiles appear in Figure 8. 

V. Electric Current Distribution 

The data have now been compiled in a form sufficient for the 
calculation of the electric current distribution from equation 23. 

One must keep in mind that the current at a particular L, cp is actually 
a function of the local partial pressures p g and p n (equation 7), but 
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we have replaced the problem of calculating them by utilizing the parameters . 
a (L) and ^ (L). The set of these parameters which is applicable to the 
radiation at the point of calculation of i must be used. Since a proton 
with an equatorial pitch angle of 30° mirrors at a magnetic latitude of 32°, 
the calculation of i for latitudes less than 32° will employ the set of parameters 
describing protons with equatorial pitch angles from 30° to 90°, and the 
calculation of i for latitudes greater than 32® will utilize the parameters 
for protons of pitch angles 0° to 30°. There will be a discontinuity in i 
as a function of latitude at 32° because as this latitude is approached from 
the equatorial side, the parameters less accurately describe the actual 
radiation. 

Two other discontinuities which arise in the distribution of i over a 
meridian plane are due to the approximation of a (L) and (L) by the simple 
mathematical expressions listed in Table III. Except for Jj (L) in the region 
L = 2.25 to 2.75 for protons with 30° to 90° pitch angles, and L = 2.75 to 
3.25 for 0° to 30° pitch angles, both a and could be well approximated 
by either straight lines or exponentials. While the approximation curves 
are continuous, their derivatives contain discontinuities between the 
segments of the curves. Hence, since i depends upon the derivatives, it 
will contain discontinuities. 

The calculations of i in both pitch angles ranges have been performed 
for the region of L from 1.5 to 10 in intervals of 0.1 and latitude from 
the equator down an L shell to approximately 1000 km above the surface of 
the earth in 2° intervals by utilizing an IBM 7094 computer at Goddard Space 
Flight Center. Actually the calculations were terminated at a latitude 



-19- 


within 2° of the latitude where the line of force attained a height of 1000 km. 
This convention was adopted so that a constant latitude increment could be used 
in the numerical integration involved in the magnetic field calculations. 

A contour plot of the current distribution on a magnetic meridian plane appears 
in Figure 9. 

The contributions of each of the terms in equation 23 to the total 
current are displayed in Figure 10 for two different L shells, one on the 
inner side of the belt where the total current is eastward, the other on 
the outer side of the belt where the total current is westward. One notices 
from an inspection of Figures 9 and 10 that the westward portion of the 
ring current, producing a southward perturbing magnetic field on the earth's 
surface, arises from the second term J 2 , which depends upon the gradiant 
of the momentum flux on the outer side of the belt. 

Since the discontinuities in i appearing in these figures originate 
in the mathematical expressions used to fit the parameters appearing in the 
equation for i, a smoothed version of the current contour plot to illustrate 
the dominant features of the current distribution appears in Figure 11. 

An integration has been performed to obtain the total eastward and 
westward currents, their algebraic sum, and the magnetic moment of the ring 
current. These values are tabulated in Table IV. 

VI. Magnetic Field Of The Ring Current. 

The magnetic field due to the ring current may be denoted by its 
components H and H respectively perpendicular and parallel to the dipole 

r Z 

axis. These components refer to a point (r, cp) specified by its polar 
coordinates (cp is latitude) or alternatively by a point (af, cp) in L space 
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where f = L cos cp. 

The ring current region is then divided into elements dS^ specified 
by (r 1 , cp') in polar coordinates or by (af',cp’) in L space, where f’ = 

L'cos cp' . Then ( Stratton, 1941), 

H = - 2 

P (acf*coscp) 

+ (2f *f 'coscp’coscp' ) Silsll I (i/F) dS 3 (25) 

‘ F 2 4 

and 

- E(k 2 ) + 2f 'coscp' (f 1 coscp' - f coscp) J (i/F)dSg 

(26) 

2 2 

Here K(k ) and E(k ) denote the complete elliptic integrals of the first and 
second kind respectively, and 



F_ 2 — f^ + ft^ „ 2f»f'cos(cp - cp’) , 

F 2 = f 2 + f' 2 + 2f.f'cos(cp + cp'). 

An element of cross sectional area of the ring current region in the 

2 A 

meridian, dS^, is expressed as dS 3 = a L ' cos^p ' dL 'dcp ' . 

The foregoing formulae appear as in Akasofu and Chapman (1961) except 
for the correction of several misprints. 

The components H p and H z were calculated separately at each point in 
the specified current region for particles in both pitch angle intervals. 


JJ (f sincp - f'sincp') [ E(k 2 ) - K(k 2 ) 
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The results were then combined to give the total field components at each 
point in the current region for all ring current particles. 

The integrations involved in calculating these field components were 
performed numerically on the IBM 7094 computer at Goddard, and the results 
appear in Figure 12. As with the current calculation, increments of 0.1 
in L and 2° (0.0349 radian) in latitude were used. 

A discontinuity in the components H and H arises when L = L f and 

P z 

2 

cp — cp’> since in this case F_ = 0, which appears in the denominator of 
the argument for each field component* The contribution to the components 
from the element of integration causing this discontinuity was taken to be 
zero each time such a discontinuity appeared. 

The coordinate system used indicates a negative H is in a southerly 

z 

direction and a negative points toward the earth. 

VII. Discussion., 

(a) Vector Field of the Ring Current. 

At the magnetic equator on the surface of the earth the magnetic 
field due to the ring current attained a value of 9 y in a southern 
direction, a value considerably smaller than the 38 y obtained by 
Akasofu , Cain and Chapman (1962) from the preliminary data from the same 
detector. A reasonable error analysis on the field strength appears 
virtually impossible to perform. Purely an estimate would place a 
confidence of about 50% (+5 y) on the field calculations. 

The distribution of field vectors from the ring current which appears 
in Figure 12 indicates a field fairly uniform on the inside of the proton 
belt and parallel to the dipole axis, as expected. In the region of the 
belt the field shows considerable curl and achieves a maximum value of 
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23 y on the equator at an L of 3.6. Outside the belt the field becomes 
parallel to the earth's field, but has the extremely small value of 2.5 y. 
Hence, the ring current field due to these protons above an energy of 100 kev 
would have little influence on the shape of the boundary of the magnetosphere. 

(b) Second Order Calculation. 

While these protons with energies above 100 kev are showing some 
magnetic influence, the magnitude is not sufficient to perturb appreciably 
the magnetic field of the magnetosphere which arises from sources internal 
to the earth. In fact it appears that the distortion of the field even at 
the altitudes of the ring current due to the compression of the field by the 
solar wind is comparable to the disturbance from the ring current ( Cahill 
and Amazeen , 1963). Hence the assumption of analyzing the data in the L, 
equatorial pitch angle coordinate system for the ring current calculations 
will produce results as accurate and physically meaningful as would an 
attempt at a second order calculation ( Beard , 1962; Akasofu , Cain and 
Chapman , 1961), in which the trapped protons are assumed to move in a 
magnetic field differing from the internal field by the field of the ring 
current itself. 

(c) Cut-off Rigidity of Solar Protons. 

In order to explain the observations of low-rigidity solar protons 

%* 

near the earth at latitudes well below those predicted by normal Stormer 
theory, investigators have included the magnetic effects of a ring current 
as well as the termination of the geomagnetic field ( Ke 1 1 ogg and Wir.ckler , 
1961; Akasofu, Lin and Van Allen , 1963; Ray, 1964). In order that a ring 
current contribute significantly to a reduction in the cut-off rigidity. 


-23- 


its magnetic moment must be some appreciable fraction of the magnetic 
moment of the earth. However, the moment due to the quiet time protons 
with energies above 100 kev is calculated to be only 0.029 that of the 
earth, a value much too small to affect solar proton trajectories. 

(d) Model Belt Parameters. 

On the basis of some very preliminary data from the proton detector 

aboard Explorer XII, the magnetic affects of a model belt that resembled 

the measured belt were calculated ( Akasofu , Cain and Chapman , 1962), 

Since then the properties of a similar belt have been applied to various 

problems of a geophysical nature (see especially the later articles by 

Akasofu; Vestine , 1963). For this work the following parameters have 

been utilized to characterize the distribution of protons: 

r Q : distance where the number density attains its maximum value, 

g^ : Gaussian parameter for the inner side of the belt: 

n = n o exp(-g^2 g 2^ where z = (r - r Q )/a, 

§2 : Gaussian parameter for the outer side of the belt, 

or. pitch angle parameter, 

n Q : maximum number density, 

E: average energy of all protons in the belt. 

In table V are listed the values adopted in the literature and a set of 

suggested values which more accurately produce the results obtained in this 

more thorough study. 

(e) Lower Energy Protons. 

One must keep clearly in mind in this discussion of the magnetic 
properties of the trapped protons that only those with energies above 

•*r. 

100 kev, the low energy cutoff for the direct detection of proton 
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particle fluxes by the detector, have been considered. However, in another 
mode of operation the detector could measure proton energy fluxes down to 
energies of some tens of kev, depending upon the steepness of the spectrum. 
Beginning at about an L of 5, large fluxes of protons were apparently 
present with energies extending down to at least 50 kev, and the ratio of 
the flux between 50 and 100 kev to that above 100 kev became increasingly 
large with range. 

While a complete saturation of the field beyond five earth radii with 
very low energy protons would significantly contribute to a quiet time ring 
current, the maximum field produced at the earth's surface by such a 
particle distribution is insufficient to account for large magnetic storms. 

For a rough estimation of this maximum field first assume that the 
particle energy density is limited by the energy density of the field on 
the magnetic equator of an L shell. Next assume that the particle intensity 
is isotropic (a = 0) and hence constant everywhere on an L shell to 
latitudes of 60°, which is essentially the surface of the earth for 
L > 5. Using equation 23 with l/hj add 1/R C evaluated at the equator, 
their maximum values, and with 

H 2 

* (L) = 2 ° 

8 tt L 6 

the current is 

i(L) = C **° cos^ep 

2 rr a 


I 
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This produces a field at the center of the earth, and approximately at the 
surface, of 


H 


z 


lxlO 4 

X 3 


Y> 


where JC. is the minimum distance of the belt from the center of the earth. 

For £= 5, H is -76 y, which is one, and probably two, factors of two too 
z 

large if compared with a field derived from any realistic distribution of 
particles exactly calculated. 

To account for a large storm of 300 y (excluding inductive effects in 
the earth) the field would have to be saturated down to at least an L of 
3.2, and probably more like 2.5 for a more reasonable particle distribution. 
The apparent necessity for such a high intensity, low altitude belt leads one 
to question whether the current equation 7 includes all the necessary terms 
to account for the ring current as it naturally occurs. One is tempted to 
suggest that electric fields may assist in the production of an increased 
westward movement of the protons already observed. It certainly would seem 
reasonable to attempt a measurement of electric fields in this region of 
space. 


It is interesting to note that the lifetime of 50 kev protons at 5 earth 
radii is only on the order of a week ( Liemohn , 1961). In order to sustain 
a quiet time belt of such particles it would be necessary to have a source 
which operates continually. 

In passing it is perhaps well to mention that the trapped electrons 
observed to date produce negligible magnetic fields. Since the magnetic 
properties of a particle depend upon the square root of its mass, electron 
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fluxes of the order of 10 electrons/cm -sec would be required to produce 
magnetic effects equivalent to those of the protons observed. This flux is 
two orders of magnitude greater than any reported in the literature. 

VIII. Conclusions. 

The protons with energies above 100 kev observed by the low energy 
proton detector flown by Davis and Williamson aboard Explorer XII are 
the first group of particles discovered in the trapped radiation region 
with intensities sufficient to produce significant magnetic effects. 
However, for applications to other geophysical problems, such as the 
distortion of field lines, the change in the Stormer cut-off rigidities, 
and the quiet time remains of a storm time ring current, the observed 
particle intensities are probably disappointingly small. 

For future studies of the proton ring current two improvements in the 
data acquisition would be desirable. The orbit of the satellite should be 
more nearly equatorial over all L values, and the detectors should be 
capable of measuring protons down to nearly one kev in energy. 

While this detector has provided information which has 
advanced our knowledge of the ring current problem, it appears that the 
complete understanding of the phenomenon will require some years of 
continued research with new particle detectors and possible complementary 
measurements of electric fields or a more complete understanding of other 


source mechaniams for the electric currents 
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APPENDIX 


Notation 

a = earth's radius, in cm. 
c = speed of light, in cm/sec. 

r = radial distance from the center of the earth to a point of interest, 
in cm. 

L = Mcllwain's coordinate, in earth radii. 

B = magnetic field vector, in gauss. 

2 

P m = magnetic pressure = B /8 tt« 
m = mass of a trapped particle in grams. 

e = magnitude of charge of a trapped particle, in esu. 
w = velocity of a particle, in cm/sec. 
w g = component of w along B. 
w R = component of w normal to B. 

E = kinetic energy of a particle, in kev. 

N = number density of particles. 

9 = pitch angle of a trapped particle, i.e. the angle between w and B, 
or its complement if larger than 90°. 

P s = pressure of the gas (trapped particles) along B. 

P n = pressure of the gas normal to B, 

cp = geomagnetic latitude. 

Note: the subscript "e" indicates that a quantity is being evaluated at 

the magnetic equator. 


TABLE I 


PROTON ENERGY SPECTRA 


30° - 90° Equatorial Pitch Angles 0° - 30° Equatorial Pitch Angles 


L 

N 0I 

E o 

E' 

E 1 

N oi 

E o 

E' 

E 1 

1.8 

POWER LAW 



POWER LAW 



2.0 

SPECTRA 



SPECTRA 



2.2 




662 (1) 

4.69xl0 7 

625 

1650 

_(3&4) 

2.5 

2.04xl0 8 

635 

1000 

__(2) 

1.46x10® 

515 

1650 

(3) 

2.8 

2.17x10® 

456 

1000 

(2) 

3 . 07x10® 

390 

900 

725 (5) 

3.0 

4.20xl0 8 

303 

900 

550 (3 > 

4.15x10® 

255 

900 

320 ^ 

3.25 

5.07xl0 8 

222 

1000 

480 

4.22x10® 

205 

900 

480 

3.5 

6.10xl0 8 

193 

880 

300 

4.64x10® 

185 

900 

226 

3.7 

4.88x10® 

196 

880 

196 

3.40x10® 

176 

OO 

— 

3.9 

3.60x10® 

195 

880 

247 

2.50x10® 

145 

900 

260 

4.6 

3.00x10® 

110 

255 

187 

1.68x10® 

105 

OO 

— 

5.0 

2.80x10® 

73 

255 

139 

1.73x10® 

61.5 

255 

90 

5.5 

1 . 62x10® 

68 

350 

148 

1.14x10® 

52 

OO 

— 

6.0 

2.05x10® 

50 

375 

275 

2.20x10® 

36 

OO 

— 

7.05 

1.35x10® 

33.5 

•» — 

OO 

6.13xl0 7 

33 

00 

— _ 


Notes: All energies in kev. 

(1) E^ = 662 kev from 470 to 1000 kev; for E > 1000 kev, power 

law with y for integral spectrum = 4.2. From 140 to 470 kev, 

no protons. 

(2) For E > 1000 kev, power law with y = 4.2. 

(3) For E > 1650 kev, power law with y = 4.2. 

(4) Below E = 255 kev, no protons. 

(5) For E > 2000 kev, power law with y = 4.2. 


TABLE II 


E c (kev) 

e c 

a 

E 0 (kev) 

N 0I 

105 

40° 

2.5 

222 

5.01x10' 

470 

30° 

2.5 

222 

5.34x10 

880 

0 

O 

2.5 

222 

4.87x10 



TABLE III 
30° to 90° E.P.A 


a 

L Region 

-Ju 

L Region 

-1.07L + 6.36 

1.5 

to 3.7 

2.95xl0' 11 exp(L/0.256) 

1.5 

to 2.4 

-0.12L + 2.84 

3.7 

to 10.0 

1. 98x10 " 7 exp(L/4.12) 

2.7 

to 3.5 




3.50xl0' 5 exp (-L/0 .810) 

3.5 

to 10.0 
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to 30° E.P.A. 



-1.62L + 6.57 

1.5 

to 2.5 

2.06x10 " 12 exp (5.18L) 

1.5 

to 2.0 

0.52L + 1.22 

2.5 

to 3.0 

4.23x10 ^exp(2.52L) 

2.0 

to 2.75 

27. 5L" 2 * 113 

3.0 

to 5.3 

5.91x10 5 exp(-1.476L) 

3.2 

to 10.0 

-0.053L + 1.16 

5.3 

to 10.0 





In the regions of 4 not included by these formulae, tabulated values from 
the curves of Figure 8 were used. 


TABLE IV 



esu/sec 

amps 

Eastward Current 

0.57xl0 15 

0.19x10 

Westward Current 

2.34xl0 15 

0.78x10 

Total Current 

1.77xl0 15 

0.59x10 


25 

Magnetic Moment 0*029M^ (M^, = 8.1x10 


3 . 

gauss -cm ) 



TABLE V 


Revised 

Value 


Value in Current 
Parameter Literature 


r o 

3.0 and 3.2 

3.5 

«1 

2.990 

1.0 

g 2 

0.419 

0.90 

a 

2.0 

2.5 

n o 

0. 6/cc 

0.46/cc 

E 

500 kev 

300 kev 
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FIGURE CAPTIONS 


Figure la. The source of the current due to the second term in equation 7. 
The eastward current between the two curved lines of force is caused by 
the crowding together of the circular orbits on the concave sides of the 
lines • 


Figure lb. The source of the current arising from the last term in equation 
7. If the particle intensity about the line of force further away from 
the earth exceeds that about the closer line, indicated by the relative 
weights of the gyration circles, a net eastward current results in the 
region between the two lines of force. 

Figure 2. The projection on the earth ! s equatorial plane of pass number 
10 on August 27 and 28, 1961. 

Figure 3. The planetary magnetic three -hour-range index Kp for the period 
during which the data used in this study were acquired. 

Figure 4. Double sloped integral proton spectrum employed in the calculation 
of the momentum integral. 

Figure 5. Sample fits of the equatorial pitch angle distributions to the 
function sin a for the regions 90° to 30° and 30° to 0°. 


Figure 6. The pitch angle parameter a as a function of energy for L values 
of 2 and 4. 

Figure 7. The values of a as a function of L for the two pitch angle ranges 
derived for the particles having the average momentum on the given L shell. 

Figure 8. The omnidirectional momentum flux (L^ for each of the two pitch 
angle ranges in units of gm/cm-sec (or ergs/cnr). , 


Figure 9. A contour plot of the total electric current in a magnetic meridian 
plane derived for protons with energies above 100 kev. The currents, 


expressed as XY signify X®10y. 
amps/km^, multiply by three. 


To convert currents from esu/cm -sec to 


Figure 10. The contributions of each of the terms in equation 23 to the total 
current for two different L shells. 


Figure 11. A smoothed version of the electric current contour plot to 
illustrate the dominant features of the electric current distribution. 

Figure 12. Vector magnetic field on a meridian plane due to the proton 
ring current. 
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